Inflation and accelerated expansion TeVeS cosmological solutions 



Luz Maria Diaz-Rivera, Lado Samushia, and Bharat Ratra 
Department of Physics, Kansas State University, 116 Card-well Hall, Manhattan, KS 66506, USA 

(Dated: January 2006 KSUPT-06/1) 

We find exact exponentially expanding and contracting de Sitter solutions of the spatially ho- 
mogeneous TeVeS cosmological equations of motion in the vacuum TeVeS model and a power law 
accelerated expanding solution in the presence of an additional ideal fluid with equation of state pa- 
rameter —5/3 < uj < — f . A preliminary stability analysis shows that the expanding vacuum solution 
is stable, while in the ideal fluid case stability depends on model parameter values. These solutions 
might provide a basis for incorporating early-time inflation or late-time accelerated expansion in 
TeVeS cosmology. 
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(N ! I. INTRODUCTION 

For more than seven decades hypothetical dark matter has been postulated to explain observed large-scale velocities 
that are larger than expected on the basis of Newton's second law of motion, Newton's inverse-square law of gravitation, 
and the assumption that cosmological structures under investigation are in gravitational equilibrium, []J, y, |3| ■ Great 
effort has been devoted to the search for direct evidence of dark matter, Q, and although there have been advances 
there is no conclusive direct evidence for dark matter. 

Observations disfavor the alternate possibility that the gravitational force falls off slower with distance than New- 
ton's inverse-square law predicts and that there is not a significant amount of dark matter (see Sees. IV. A. 1 and 
IV.B.13 of Ref. @). 

Another explanation — modified Newtonian dynamics (MOND) — has been proposed by Milgrom Q. In MOND, 
Newton's second law, a = — V^at, is modified to 
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Here a is acceleration, do is an acceleration scale, $jv is the Newtonian potential, and n(x) is a function which 
satisfies /j,(x) » 1 when x ^> 1, and (i(x) ~ x when x -C 1. From large-scale galactic data it has been estimated 
that ao w 1 x 1CU 8 cms -2 , so in the solar system where accelerations are large compared to ao the usual Newtonian 
+3 , law is recovered. Equation and this choice of fJ,(x) was proposed to explain the observed flat rotation curves of 
disk galaxies, without need for dark matter. MOND also explains the Tully-Fisher relation, the observed correlation 
between the infrared luminosity of a disk galaxy and the rotation velocity in the flat region of the rotation curve. See 
Refs. for other tests and predictions of MOND. 

In spite of its success, MOND is not a consistent theory and has some apparent problems. MOND requires dark 
matter in galaxy clusters, 0, and if this is massive neutrinos the needed neutrino mass is on the edge of being ruled 
out 0- There are indications that observations require that the value of ao depend on the object being studied [Tof . 
Other possible problems are considered in Refs. [IlJ. A major drawback of MOND was the lack of a relativistic 
generalization needed for a consistent MOND cosmological model. 

Bekenstein 0] nas recently proposed a new covariant field theory which has MOND characteristics in the weak 
acceleration limit and provides a setting for constructing consistent cosmological models. This model uses three 
dynamical gravitational fields — a tensor, a vector, and a scalar field — leading to the acronym TeVeS. Several authors 
have studied different aspects of the TeVeS model, including large-scale structure formation [13( , gravitational lensing 
[T^ . and the strong gravity regime of the model Hao and Akhoury have studied a cosmological model and 
its dynamics in the context of TeVeS , but no exact solution has been found so far. 

Whether TeVeS or some generalization is able to provide an accurate model of cosmology is an issue of some 
current interest. The standard dark energy and dark matter dominated general relativistic cosmological model does a 
remarkable job on large scales, however, it appears to have trouble fitting some smaller scale observations TeVeS 
cosmology on the other hand is only now attracting preliminary attention. In the standard cosmology the Universe 
goes through two periods of accelerated expansion, at early times during inflation and during the current epoch. 
Inflation is used to make the Universe homogeneous and to generate the quantum-mechanical fluctuations responsible 
for observed large-scale structures ^3 wnue supernova and other data indicate that the expansion is speeding up now 
It is of interest to understand whether TeVeS cosmology has an accelerated solution that might play a similar 
role to that in the standard model. A complete and consistent TeVeS cosmological model would be very useful, at 
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the very least providing stimulus for the development of precision tests that can be used to distinguish the TeVeS 
predictions from those of standard cosmology, along the lines of observational tests that were developed to distinguish 
open inflation from the standard flat model |18j . 

The paper is organized as follows. In Sec. [H] we present the general TeVeS equations of motion and in Sec. 
IIIII we specialize these equations to the spatially homogeneous and isotropic Friedmann-Lemaitre- Robertson- Walker 
spacetime model. In Sec. IIVI we present exact solutions of these equations and discuss cosmological implications. 
We examine the stability of these solutions in Sec. [V] Finally, in Sec. IVII we summarize our results and present 
conclusions. 



II. TEVES EQUATIONS OF MOTION 

TeVeS assumes three dynamical gravitational fields: a tensor metric g^, a unit norm time-like four- vector U a , and 
a scalar <j>. There is also a nondynamical scalar field a. The three dynamical fields are connected through the physical 
metric tensor 



g al3 = e~ 2 *g al3 - 2U a Up sinh(20) = e" 2 *^ + U a U p ) - e 2 ^U a U p , 



(2) 



where g a p is the Einstein metric tensor. Our conventions are those of Bekenstein [l2j with metric signature 
diag(-l, 1,1,1). 

Following Bekenstein the action of the system is written as 



S Sg -\- Ss ~\~ Sy "I - Sm-i 



(3) 



where the actions for the tensor field S g , the two scalar fields S s , the vector field S v , and ordinary matter (as opposed 
to Bekenstein's <f>, a, and U a fields) S m , are, 
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Here g and R a p are the determinant and Ricci tensor constructed from the Einstein metric tensor g a p, the tensor 
h a @ = g al3 — U a lJP has been used instead of g a ^ in order to avoid acausal dynamical scalar field propagation [l9|. 
and f a represents ordinary matter, including possibly a scalar field or a cosmological constant that drives accelerated 
expansion or inflation. Also, F is a free dimensionless function, K and k are positive dimensionless coupling constants, 
I is a constant scale of length, A is a spacetime dependent Lagrange multiplier, and G is the bare gravitational constant. 
Square brackets in Eq. © denote antisymmetrization, Ay^B^ = A^B^ — A^B^, in what follows round brackets denote 
symmetrization, A^B^ = A^B V + A V B^, and the covariant derivative denoted by "|" in Eq. Q is defined in terms 
of the physical metric tensor. 

The equations of motion are obtained by varying the action given in Eq. @ with respect to U a , cr, <j), g a p, and A. 
These equations are, 

K{U^. fi + U a U 7 U^ ;l3 ) + %nGo 2 \UP<\> 3 g a "t<\> a + U a {U^4>^f] = 8%G(1 - er^g^U^f^ + U a U fi XPf^p], (8) 



- -^ 2 F'(/i) = kl 2 h al3 4> l0t ^, 



Hkfh^^^^h^^y = kG[g aS3 + (1 + e-^)U a U fi ]f aP , 



(9) 
(10) 
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G a(3 = 8vG[T aP + (1 - e~ 4 *)U' t T^ a U 0) + r Q/3 ] + B a0 , (11) 



g^U M U v = -l. (12) 

Here T^ v is the stress-energy tensor of ordinary matter. Equation (JHJ is the vector field equation of motion. The 
prime in Eq. © denotes a derivative with respect to p, = kGa 2 . Equation ® determines p and the nondynamical 
scalar field a in terms of g a /3i U a , and (f>. Equations (|10f) and Ijlll) are the equations of motion for the dynamical 
scalar field <p an d the metric tensor respectively. The Einstein tensor G a p = R a p — Rg a p/2, where R is the Ricci 
scalar constructed from the metric tensor g a p. We call the right hand side of Eq. lfTT|) 8ttGT^. 
The tensors r a p and Q a p in Eq. i|ll|) are defined as 
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^F(p)g a0 , (13) 



©a/3 = k( 9 ^U m U m - -g^g^U^U^g^y (14) 

Assuming that ordinary matter is an ideal fluid, the matter stress-energy tensor in the physical coordinate system is 
T a = pu a ii/3 +p(g a f3 + UasUp), with p the energy density p the pressure, and u a the fluid four- velocity. An equation 
of state p(p) must be specified to fully determine this stress-energy tensor. The Bianchi identity T^i ^ = 0, together 
with Eqs. ||S1)- (|1U|) and (|12|l . implies stress-energy conservation for ordinary matter, which may be used as an equation 
of motion for ordinary matter, instead of using part of the tensor equation of motion, Eq. (|lll) . Equation (|12fl . which 
is a constraint equation derived by varying the action with respect to the Lagrange multiplier A, forces U a to be 
time-like with unit norm. 

III. SPATIALLY HOMOGENEOUS AND ISOTROPIC EQUATIONS OF MOTION 

Observations indicate that on large scales the spatial distribution of matter and radiation are close to isotropic 
in the mean Q ■ This fact motivates consideration of the spatially homogeneous and isotropic Friedmann-Lemaitre- 
Robertson- Walker spacetime model. In general relativistic cosmology observations also show that space is close to 
flat (see, e.g., 0). Although it is not yet clear if observations also favor vanishing spatial curvature in Bekenstein's 
model, for simplicity we will consider the spatially-flat Friedmann-Lemaitre-Robertson- Walker metric, 

dS 2 = -dP + a 2 [dr 2 +r 2 (d6 2 + sin 2 d dtp 2 )}, (15) 

where dt = e^dt and a — e~^a relate the time and scale factor in the physical and Einstein coordinates. The scalar 
field cj> depends on time and since there cannot be a preferred spatial direction the vector field U 11 must point in the 
time direction, Ui 1 = 5% [lj. 

In the spatially homogeneous and isotropic model, the tensor field equations of motion obtained from Eq. I|ll|) are 

a 2 -— PS ' + Vk^ + 3W^< (16) 

2± + h - 2 =-8,Gpe^-^ + %F{»), (17) 

where an overdot denotes a time derivative. These equations are the analog of the Friedmann equations in general 
relativistic cosmology. The dynamical scalar field equation of motion obtained from Eq. (|10|l is 

(i4>+ (3// + /i)0+^e- 2 *G5 + 3p)=O. (18) 
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Additionally, stress-energy conservation obtained from the Bianchi identity, T^L = 0, and Eqs. l(5|)- (|10|) and l|12() 
leads to the equation of conservation of stress-energy of ordinary matter, 

P = 3(<j>-^)(p + p). (19) 

Consistent with the case of general relativistic cosmology, it may be shown that only three of the four equations 
H16f) - I|19|l are independent. Note that C/ M = 5^ causes U^'^ to vanish, so Eqs. (JSJl and l|12l) are identically satisfied. 
On the other hand, given an expression for the free function F{p) (which must be externally specified since there is 
no theory for its determination), Eq. determines p and a in terms of <f), U a , and g a p- Equations ©, @, and l|12|) 
may be ignored in the following analysis of the spatially homogeneous and isotropic model. 

In the following section we look for analytical solutions of the set of field equations l|16fl and (|18|l along with (|19|l . 

IV. SOLUTIONS 
A. Vacuum case 

The vacuum case with p = = p is a simple case to solve analytically. The Friedmann and scalar field equations 
become 

a 2 = Tk^ + ^ F ^ (20) 

n4> + + frjjt = 0. (21) 

Assuming a time independent p, we find two exact solutions. 

1/2 

Our first solution corresponds to a = ciq and <j)(t) — at + <j)Q, where a = [— pF / (4kl 2 )\ . Since the scale factor is 
constant, this solution is not very interesting. 

Our second solution is the analog of the general relativistic de Sitter solution 

a(t) = a e ±Ho \ (22) 
<t>{t) = 0o, (23) 

1 /2 

where a and 4>o are constants and H = [2ir p 2 F '/ '(3k 2 l 2 )] . This solution also depends on the form of the free 
function F(p), which is required to be positive in this case. This solution corresponds to the "slow roll" approximation 
solution discussed in Ref. [l^. It is interesting that the TeVeS model has an accelerated expansion solution even in 
the absence of a cosmological constant or a dark energy scalar field. 

B. Barotropic fluid case 

We assume that the ordinary matter fluid has a barotropic equation of state p = uip,so the stress-energy conservation 
equation ijT§|) has the solution 

e 3(l+w)0 

In order to find an analytical solution of the field equations, we make the ansatz 

[ o ( t )]3(l+o0 = e (l+3 W )0« (25) 
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so pit) — exp(20). The Friedmann and scalar field equations, (|16fl and (|18(l . then become 



1 (1 + 3c^) 2 
3 (1+w) 2 
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fxcj) + -r—^ 2 +^+—{l + iuj) = 0. 
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From Eq. Il2(ill we note that if p, is time independent, <j> is time independent and so 
simplifies Eq. I|27l) . resulting in an analytical expression for the scalar field, 



(27) 

0. Taking <j> — considerably 



4>(t) = ±at + (f> , 



(28) 
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where <f>o is a constant of integration and now a — [— fcG(l + w)/(2/i)] . Using this solution for cj)(t) in the Friedmann 
equation l(2f)j) we get a relation between the constants u> and p and the free function F{p), 
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/x 2 F( M ) 



fcG(l + 3u;) ; 
6^(1 + cj) 



4ttG(1 -uj). 



(29) 



It is important to note that this equation does not give us a function -F(/i), but instead given an F{p) it relates fj, 
and uj, or given an F(fi) and an uj it fixes the value of /i. 

From Eqs. (|25() and 128(1 we see that our solution is a de Sitter solution with scale factor 



ait) = age 



±H t 



(30) 



where 



Ho = \ — 



fcG(l + 3w) 2 



18/i(l + w) 



(31) 
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(32) 
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Recalling the definition = fcGcr 2 , we get Ho = (1 + [— I8cr 2 (l +w)] . Since i? must be real, our solution 
must have uj < — 1. This is not necessarily a problem, since we are interested in a solution that may be applicable 
during inflation or the current accelerated expansion epoch. It is interesting that Ea. H3U|) is a de Sitter solution, 
independent of the value of uj, provided w < — 1. 

Using this solution in the field equations, we find the following useful identities, 



8ttG 



1 + W- + Al^G FM 



(33) 



3H a 
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3lu) = 0. 



(34) 



The explicit contribution of the scalar field to the energy density in this model is not easily seen because ordinary 
matter and the scalar field are entangled, not only in the metric, Eq. (J2J, but also in the right hand side of Eq. I|II|) . 
f 00 = pe 2<#> is the energy density of ordinary matter so f G ff - Too = -pe 2 ^ + pe~ 20 + ^ 2 /(fcG) + /i 2 T/(4fc 2 Z 2 Gjcould 
be associated with the contribution of the scalar field to the energy density, as in Brans-Dicke-like models 20] , but 
note that in Tqq — Too the stress-energy tensor of ordinary matter also contributes. Of course, Tq^ is the right hand 
side of the Friedmann equation which we rewrite as 



8ttG 



-Poff, 



(35) 
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and similarly the spatial components Tff on the right hand side of Eq. 1|17|) ma y be rewritten as 

[Pes + 3p c ff] , (36) 



a 4ttG 



a 

where 

Peg and p c g obey the energy conservation equation p e g — — 3(/2 e g + p c s)a/ 'a. It may be shown that when p e g and 
p c g, Eqs. I|37|l and are used in this spatially homogeneous conservation equation, we recover the spatially 

homogeneous part of the scalar field equation of motion, Eq. H18f) . See, e.g., Ref. ,21] for the general relativistic 
cosmology analog of this result. 

Note that in Ref. 0] p m = pe~ 2 ^ is defined to be the energy density of ordinary matter, and p^ = p4> 2 /(kG) + 
p 2 F/ (Ak 2 l 2 G) the effective energy density of the scalar field 4>. However, we have shown here that the energy density of 
ordinary matter is given by the time-time component of the stress-energy tensor of ordinary matter p m = Too = pe 2 ^. 
The definition of p^ in Ref. corresponds to a piece of T§ff — Too-that may be associated with the energy density 
of the scalar field. Thus, the energy densities p m and p^ defined in [lfj do not correspond to the energy densities of 
a fluid in the sense that they do not obey an energy conservation equation obtained from the Bianchi identity, which 
reduces to an equation of motion for the scalar field. These definitions are however, a piece of the effective energy 
density defined above. 

Our solution, Eqs. 1|28|1 -H30 [1 . implies that p c e = —k(l + 3w) 2 /[487r^(l + ui)] is positive for to < —1. 
We will examine the stability of our solution in more detail in the following section. 

C. Bekenstein's F(p) function 

The exact solutions that we have obtained are valid for any F(p) function, as long as the function satisfies the 
relation given in Eq. (|29|l for the solution corresponding to a barotropic equation of state, Sec. IIVBI and F(p) > 
for the solution corresponding to the vacuum case, Sec. IIV Al As pointed out by Bekenstein, [12J, there is no theory 
for the function F(p), and based on this freedom and other considerations he picks as an example, 

F(p) = ^-(4 + 2p Ap 2 + p 3 ) + JL ln[(l - p,) 2 ]. (39) 

The cosmologically relevant part of this function lies in the interval 2 < p < oo. Using this expression, from Eq. 
(12911 we confirm that our solution corresponding to a barotropic perfect fluid is limited to to < — 1. We also find that 
-F(ju) > for any value of p in the interval relevant to cosmology, so this form of F(p) can also be used in the solution 
for the vacuum case. 

D. Density and deceleration parameters 

We rewrite the Friedmann equation, Eq. (|16|l . as 

1 = fi m + fy, (40) 
where the density parameters of ordinary matter and the scalar field are, 

= ^pe 2 *, (41) 



3H 



3H 2 



pe^ + pe- 2 * + ^ + 7 £_F(p) , (42) 



kG^ ' Ak 2 l 2 G 
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where H = a/a = H + <fi e ^, H = (da/dt)/a. As in general relativistic cosmology, we use Eqs. <|16ll an d l|17|l to 
determine the deceleration parameter, q = —a/aH 2 , finding 



47rG _ 87rG 



3ff 2 



3iJ 2 



2^ ; 2 
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(43) 



For our barotropic perfect fluid solution, Eqs. I|28[) - l|32() . Eq. (|43() gives q = — 1 independent of the form of the -F(/Lt) 
function and the value of to. The same result also holds for our solution in the vacuum case (Sec. IIV Af) . Of course, 
both results follow directly from Eqs. (|2"2"|l and (fTTTTfl and the definition of q. 



E. Physical variables 

Until this point we have been working in the Einstein metric, i.e., in the metric g M „ whose dynamics is governed 
by the Einstein-Hilbert action. We will call this frame the "Einstein frame". It is important to write the solutions 
found above in terms of the physical variables, i.e., in the frame of the metric g^ v known as "physical frame" or 
what Skordis et al. called the "matter frame". These metrics are related by Eq. J2J, so for the spatially- flat 
Friedmann-Lemaitre-Robertson- Walker metric of Eq. I|15f) we have a = e~^a and dt = e^dt. In the "physical frame" 
we define the deceleration parameter in terms of physical variables as 



Q = 



a(i)a(t) 



(44) 



where now an overdot denotes differentiation with respect to physical time t. 
• Vacuum case. 

The solution of the vacuum case in the Einstein frame, Eqs. and l|23|) . corresponds to the physical frame 
solution 



a(t) = doe 



±H i 



(45) 



where do — aoe ^° and Ho — Hoe ^° . The vacuum solution in the physical variables is also a de Sitter solution 
with q = — 1, as can be seen from Eq. (|44|l . 

• Barotropic ideal fluid case. 

We have found that the solution for a barotropic perfect fluid in the Einstein frame is a(t) — aje °' and 
(p(t) = ±at + <fio, where H and a are given by Eqs. I|31|) and l|28|) . Rewriting this solution in terms of the 
physical variables we get 



~U\ ~ 7H a /a-l - r-2/3(l+oj) 

a(t) = oq t ' = &o t j 



(46) 



(j)(t) = \n{at) 



(47) 



where d = a e"^ - ff °/ a (a) (Ho/a_1) = [-kG(l + uj)/(2fi)]- 1 / 3( - 1+u \ Note that the contracting solution in the 
Einstein frame does not have physical meaning in the physical variables. As we can see, the de Sitter solution 
in the Einstein frame for a barotropic ideal fluid is transformed to a power law solution in the physical frame. 
Since lu < — 1 the solution is expanding. 

Using Eq. (|44|l it is straightforward to check that the deceleration parameter in the physical frame for the 
solution (|46|l is q = —(5 + 3u>)/2, then q < implies ui > —5/3. Since uj < — 1 this solution expands in an 
accelerated way for lu in the interval of values —5/3 < uj < —1. Such solutions have been discussed in the 
context of power law inflation |2J|. 
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V. STABILITY OF THE SOLUTIONS 



In Sec. II VI we have found exact solutions for the vacuum case, Eqs. (1221) and l|23|) ■ and for the case with a barotropic 
ideal fluid, Eqs. I|28|) and i|3U|) . The purpose of this section is to see if these solutions are stable with respect to small 
perturbations. 

A complete analysis will require solving the TeVeS cosmological linear perturbation equations (see Refs. [l3L 
for small departures from the spatially homogeneous and isotropic background solutions that we have derived in the 
previous section. In our preliminary exploration here we will focus only on very large scale perturbations, using 
instead the method developed in Ref. |23j . 

We first consider the barotropic ideal fluid case. Following the method of Ref. [23J we make the change of variables, 

a(t) = a e (t)v(t), 

<t>(t)=4>e{t)+m, (48) 

p(t) =p e (t)+T)(t). 

Here v(t) — 1, -0(t) , and rj{t) are small perturbations about the spatially homogeneous solution, denoted here as a e , 
(f> e , and p ei and given in Eqs. (|3"U|) . and (PH . Defining two "speeds" s — ip and r = v (note that p obeys a first 
order differential equation and so we do not need a "speed" for it), the equations of motion become, 



\3H +3~ + ~ )s e (e ^ (p e + r)) - p e ) , 

\ v p) 2p 

(49) 

(1 + 3o>)e- 2 ^ (e-^{p e +r))-p e )v- ^-v {T{a + s) - T{a)) , 





s, 




s = 


3a 
vp 


—a — 
A* 


V = 


r, 




r = 


-2H Q r - 


4irG 
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r) = 


-3(1- 


S )(l 



)p e - d^o + - — a— SJT), 

where T = p<j) 2 /(kG) + p 2 F(p)/(Al 2 k 2 G). 

The vector field U a contributes to the initial set of equations (|%|l- (|12[l . but it does not contribute to the spatially 
homogeneous cosmological equations (1 1 01) l|19|) . because there is no preferred spatial direction on cosmological scales 
and so U a — 5f . To simplify the analysis we do not consider a perturbation of U a here. 

p((j)) is a function of the scalar field's time derivative and is time independent for our background solutions p — p{a). 
When we perturb the scalar field, p = p(a + s) becomes time dependent. Since we are interested in linear perturbation 
theory, we expand it in a Taylor series with respect to the small parameter s, and we keep only terms up to first order 
p(a + s) = p(a) + s dp/ da. 

The critical point of Eqs. I|49l) is (ipo, So, t>o, ro, 770) = (0, 0, 0, 0) , where v is an arbitrary constant which corre- 
sponds to the freedom in rescaling ao- This point of phase space is where all time derivatives in Eqs. (|49f) vanish. 
Perturbing Eqs. (|49|l around the critical point, {ip, s, v, r, rf) = (jpo + ipi, sq + s±, vq + v\, ro + ri, 770 + ?/i), the linear 
perturbation equations are, 



01 


= si, 






Si 






Po 




= n, 
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3 ( 








= 3(1 + 



3a kG 
— ri-— ( 
vp 2p 



— — v Sl - 2H n - — (1 + 3Lj)e 2 ^vrn, 
6 da 3 



(50) 



< e si - 3 — + _^^ e ri - 3(H - a)(l + w)r)i, 
v 

where po — p{a). 

In order to compute the eigenvalues of the matrix defined by these equations, we need to specify the function F(p) 
and the values of the constants of the TeVeS model. As an example, we consider Bekenstein's F(p) function, 
Eq. I|39|) . for some values of p in the cosmological regime 2 < p < 00. We choose two set of arbitrary values for the 
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parameters of the model, k = I = G = 1 and k = 0.1, I = 10, G = 1. Also, we consider three illustrative values for 
a; = —1.02, —1.05, and —1.08. Numerical computation shows that the expanding solution, a(t) = aoe Hot , is stable for 
small values of [i (/i < 19 when to = —1.02, /i < 13 when w = —1.05, and /i < 11 when w = —1.08) for k = I = G = 1. 
For the second set of parameters k — 0.1, I — 10, G = 1 we find that the expanding solution is stable for small values 
of /i (p < 11) for all three values of ui that we consider. The contracting solution, a(t) — aoe~ Hot , on the other hand, 
is not stable for any /i in the cosmological allowed range, independent of the values of the other model parameters 
we consider. In the limit of small deviation from u> = — 1, analytical computations show that the expanding solution 
a(t) = aoe H ° * is stable. Clearly, the stability of the solution depends sensitively on the form of F(ji) and the values 
of k, I, G, and fi. We note that even though co < — 1, this does not necessarily imply instability. 

The expanding solution for the vacuum case, a(t) = a e Hot with (j)(t) — <f> , is stable with linear perturbation 
eigenvalues Ai = —2H , X 2 = —3H . The contracting solution, a(t) = aoe~ H ° is unstable with eigenvalues Ai = 2H Q , 
X-2 = 3H . These results hold for any F(fi) > function. 

We have calculated analytically the eigenvalues of small perturbations in the limit of general relativity. We find 
that for both the vacuum and barotropic equation of state cases, the expanding solutions a(t) — aoe Hot are stable, 
as well known. 



VI. CONCLUSIONS 

We have found two homogeneous and isotropic cosmological solutions of the TeVeS model, a relativistic general- 
ization of MOND. In the vacuum case, the analytical solutions in both, Einstein and physical frames, are de Sitter 
solutions with q = — 1, and the scalar field is a constant 4>(t) = 4>q. For a barotropic ideal fluid we find a de Sitter 
solution in the Einstein frame, which is transformed to a power law expanding solution in the physical frame. Ac- 
celerated expansion for this solution occurs for —5/3 < to < —1. The scalar field in this case evolves logarithmically 
with the physical time, (f>(t) = ln(crf). 

The free function F(/i) plays an important role in the TeVeS model. However, we have found that for vacuum case 
the only restriction is F(/i) > 0, while for the barotropic perfect fluid case our solutions are valid for a family of F(fj.) 
functions that obey the u) dependent relationship given in Eq. I|29|l . 

In the vacuum case the expanding solution a{t) = a^e Ht is stable for any F(fi) > function. In the more general 
case of a barotropic ideal fluid, we have studied stability for a few values of the parameters and constants of the 
theory and have assumed Bekenstein's F(fi) function. Typically, the expanding solution, a = a,Qe Ht in Einstein frame, 
is stable for small values of [A, and becomes unstable for large values of fi. Since the solutions in Einstein frame 
and physical frame are related through transformation of variables, the corresponding solution in the physical frame 
a(t) = aot H °/ a ~ 1 has the same stability properties. 

The accelerated expansion TeVeS cosmological solutions we have found might prove helpful when trying to incor- 
porate inflation or the present epoch of accelerated expansion in the TeVeS scenario. 
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